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Abstract 

A non-autonomous flow system is introduced with an attractor of Plykin type that 
may serve as a base for elaboration of real systems and devices demonstrating 
the structurally stable chaotic dynamics. The starting point is a map on a two- 
dimensional sphere, consisting of four stages of continuous geometrically evident 
transformations. The computations indicate that in a certain parameter range the 
map has a uniformly hyperbolic attractor. It may be represented on a plane by means 
of a stereographic projection. Accounting structural stability, a modification of the 
model is undertaken to obtain a set of two non-autonomous differential equations 
of the first order with smooth coefficients. As follows from computations, it has the 
Plykin type attractor in the Poincare cross-section. 
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In mathematical theory of dynamical systems a class of uniformly hyperbolic 
chaotic attractors is known [112,3,4, 5,6,7]. In such an attractor all orbits are 
of the saddle type, and their stable and unstable manifolds do not touch each 
other, but can only intersect transversally. These attractors manifest strong 
stochastic properties and allow detailed mathematical analysis. They are struc- 
turally stable; that means insensitivity of the structure of the attractors in 
respect to variation of functions and parameters in the dynamical equations. 
Although main concepts of the relevant mathematical theory were advanced 
40 years ago, until recently, these attractors are considered rather as a refined 
image of chaos, than as adequate models for real- world systems. In textbooks 
and reviews examples of the uniformly hyperbolic attractors traditionally are 
represented by mathematical constructions, like Plykin attractor and Smale 
- Williams solenoid. These examples relate to discrete-time systems, iterated 
maps. In particular, Plykin attractor takes place in some special map on a 
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sphere with four holes, or on a plane in a bounded domain with three holes 
[6]. In applications, physics and technology people deal usually with systems 
operating in continuous time, called flows in mathematical literature. In such 
systems Plykin type attractors could occur in the context of description based 
on the Poincare map [7f8f9] . 

The present work is aimed at explicit construction of a non-autonomous flow 
system with Plykin type attractor, which could provide a basis for develop- 
ment of real systems and devices, demonstrating structurally stable chaotic 
dynamics. The starting point is consideration of a motion on a two-dimensional 
sphere composed of periodically repeated stages of continuous geometrically 
evident transformations. 

Let us consider a sphere of unit radius. A point on the sphere can be specified 
in angular coordinates (6, ip), or in Cartesian coordinates 

x = cos^sin^, y = sin sin 6* , z = cos9, (1) 

which satisfy the relation x 2 + y 2 + z 2 = 1. As proven by Plykin, a map on 
a sphere can possess hyperbolic attractor only in the presence of at least four 
holes, the areas not visited by trajectories belonging to the attractor. In our 
construction this role will be played by neighborhoods of four points A, B, C, 
D with coordinates (x,y,z) = (±l/>/2, 0, ±l/\/2). 

Let us consider a sequence of four successive continuous transformations, each 
of which is of unit time duration. 

I. Flow down along circles of latitude, that is motion of the representative 
points on the sphere away from the meridians NABS and NDCS (N and S 
designate the north and the south poles) towards the meridians equally distant 
from the arcs AB and CD. In Cartesian coordinates it is governed by equations 

x = —exy 2 , y = ex 2 y, z = 0, (2) 

where e is a parameter. 

II. Differential rotation around z-axis with angular velocity depending on 
z linearly, in such way that the points B and C do not move, while A and D 
exchange their location; it corresponds to equations 

x = 7i(z/V2 + l/2)y, y = -7i(z/V2 + l/2)x, z = 0. (3) 

III. Flow down to the equator, that is motion of representative points 
along circles centered on the x-axis on the sphere from the great circle ABCD, 
towards the equator: 

x = 0, y = eyz 2 , z = —ey 2 z. (4) 
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IV. Differential rotation around x-axis with angular velocity depending on 
x linearly, in such way that representative points in the plane, orthogonal to 
the x-axis and containing the point C, do not move, while those in the plane 
containing the point B undergo a turn by 180°: 

x = 0, y = -7i(x/V2 + l/2)z, z = tc{x/V2 + l/2)y. (5) 



Note a symmetry of the procedure: the first and the second pairs of the trans- 
formations are identical, up to exchange of the variables x and z. 

The Poincare map describing transformation of a state vector x n = (x n , y n , z n ) 
on a period T = 4 may be obtained explicitly. From successive solution of the 
differential equations (JSl)-© with account of the mentioned symmetry, one 
can represent the resulting state vector x n+1 as 



x„+i = f+(f_(x n )), f ± (x) 



±z 



v / i 2 +j 2 (ye^+y 2 )/ 2 cos f (zV2+l)±xe-^ x2 +y 2 ^ 2 sin § 

^y x 2 e -E(x 2 +y 2 ) +y 2 e e(x 2 +y 2 ) 

^x 2 +y 2 (ye< x2 +y 2 y 2 sin f (zV2+l)Txe-^ x2 +y 2 ^ 2 cos f (W2+1)) 

^y x 2 e - £ {x 2 +y 2 ) +y 2 e e(x 2 +y 2 ) 

(6) 

The relations (jBJ) determine the map on the sphere x n+ i = T(x n ). Note that C 
is a fixed point of the map T, while A, B and D compose an unstable periodic 
orbit of period 3: A — » D — > B — > A. The map T is invertible. The inverse 
map appears as a result of the same transformations in backward order, with 
reversed directions of the rotations. 

Figure 1 shows attractor of the map T at e=0.77. Observe specific fractal-like 
transversal structure of the attractor: the object looks like composed of strips, 
each of which contains narrower strips of the next level etc. 

Description of the dynamics can be reformulated to represent instantaneous 
states of the system on a plane. The variable change is 

W = X + iY= X - Z + iy f , (7) 
x + z + Vl W 

which corresponds to a stereographic projection, with selection of the pro- 
jection point at C(— l/y/2, 0, —l/y/2). This point does not belong to the 
attractor (it is in the "hole"), so the image of the attractor on the plane is 
located in a bounded domain. Portrait of the attractor for the Poincare map 
in this representation is shown in Fig. 2a. 

I argue that the attractor relates to the uniformly hyperbolic class, i.e. it is an 
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Figure 1. Attractor of the map (J6j) at £=0.77 on the unit square 
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Figure 2. (a) Portrait of the attractor of the map (J6j) at e=0.77 presented on a plane 
by means of the stereographic projection and (b) disposition of stable (black) 
and unstable (white) manifolds in a domain D (gray background) that contains the 
attractor 

attractor of Plykin type. To support this assertion let us turn to Fig. 2b illus- 
trating location of stable and unstable manifolds in a domain D that contains 
the attractor. The stable and unstable manifolds are shown, respectively, 
by black an white curves, drawn on a gray background of the domain D. As 
seen, the unstable manifolds follow along filaments of the attractor, while the 



To draw the stable and unstable manifolds with a computer, we do the following. 
First, for a given point on the attractor x we obtain its image by iterations of 
the Poincare map x = Tp (x), and its pre-image by iterations of the inverse map 
x = Tp N (x), where N is some empirically chosen integer. Then, with random initial 
conditions y in a small neighborhood of x, we get a set of points y = Tp (y) by 
iterations of the Poincare map, which mark the unstable manifold. In a similar way, 
starting with random initial conditions y in a small neighborhood of x, we draw the 
stable manifold with a set of points y = Tp N (y). The accuracy the manifolds are 
depicted grows fast with N . Actually, N = 6 is enough to get so small errors that 
they are indistinguishable in the plot. 
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stable ones are transversal to the filaments. Mutual disposition of the stable 
and unstable manifolds certainly excludes possibility of tangencies, at least in 
the domain D. 

An alternative approach to verification of the hyperbolicity can be based on 
the cone criterion known from the mathematical literature [1,2,3,4 ,5,9J. Such 
calculations were performed on a base of methodology developed in Ref.|10|. 
and the hyperbolicity was confirmed!!] 

Lyapunov exponents were computed for the map T by means of the proce- 
dure based on joint iterations of the map together with a collection of two sets 
of linearized equations for perturbation vectors. At each iteration step, these 
two vectors are orthogonalized, and normalized to a fixed constant. Lyapunov 
exponents are obtained as slopes of the straight lines approximating the accu- 
mulating sums of logarithms of the norm ratios for the vectors in dependence of 
the number of iterations [TT]. In particular, at e=0.77 the Lyapunov exponents 
are Ai = 0.9587 and A 2 = —1.1406. Then, the estimate of the attractor di- 
mension from the Kaplan - Yorke formula yields D KY ~ 1 + A 1( /|A 2 | = 0.8405. 

One can accept an interpretation that an instant speed of a representative 
point on the sphere is determined by combination of two vector fields, which 
are switched on and off, turn by turn. One corresponds to dynamics during 
the stages of flow down, and another to the stages of differential rotation. 
Let us introduce angles a and (3 determining directions of the fields, and 
coefficients p and q responsible for switching them on and off. We set a = 
(3 = | (_i)[*/2] _ |, p = 1 + q = |(-1)[*/ 2 1 - l(-l)[*/2+i/2]_ (Here 

[r] designates the integer part of r.) Setting K = n/\/2, we write down the 
Eqs. ©-© in a compact form: 

x = —pey 2 (x cos a + z sin a) cos a + qKy(—x sin (3 + z cos (3 + cos (3, 
y = peyix cos a + z sin a) 2 — qK(x cos (3 + z sin (3){—x sin (3 + z cos (3 + ^=), 
z = ~pey 2 {x cos a + z sin a) sin a + qKy(—x sin (3 + z cos (3 + A=) sin (3. 

(8) 

As a next step, let us construct a version of the non-autonomous model con- 
taining only smooth functions. To do this, in Eqs. ([Sj) we simply set 

a = | cos ft, (3=1 sin f t, p = 1, q = sin ft. (9) 

Now, two vector fields, responsible for the flow down and for the differen- 
tial rotation vary in time continuously and smoothly, undergoing rotations in 

2 These results will be published elsewhere, because require more volume than ap- 
propriate for the short communication. 
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space in such way that the original configuration is repeated with the period 
T = 4. The second field oscillates, reversing the direction twice on a period. 
The time variations of one and other fields are shifted in phase relatively by 
a quarter of period. So, the extremal values of the fields are achieved turn by 
turn. The action of the fields at the extrema are mostly significant for mo- 
tion of the representative points, hence, it is reasonable to specify the values 
at extrema, like in the original model. Because of structural stability of the 
hyperbolic attractor, one can hope that its nature remains the same after the 
modification, at least in a properly selected range of the parameters. As seen 
from computations, it is so, e.g. at £=0.72 and K=1.9. 

In contrast to the previous version of the model, the Poincare map can not 
be expressed analytically. However, it may be easily realized by a computer 
program integrating the equations with a finite-difference method on a time 
period T = 4. 

We can exclude a redundant variable in the equations by means of the variable 
change (J7|), and rewrite them in terms of W = X + iY . After separation of 
the real and imaginary parts, we obtain 



X = -2£F 2 fi!(X, Y, t) (cos(| cos ft) - X sin(| cos ft) 
+KYn 2 (X, Y, t) (cos(f sin ft) - X sin(| sin ft)) sin f , 

(10) 

Y = eYQ^X, Y, t) (2X cos(f cos ft) + (1 - X 2 + Y 2 ) sin(f cos ft)) 
-KQ 2 (X, Y, t) (X cos(f sin ft) + f (1 - X 2 + Y 2 ) sin(f sin ft)) sin f , 

where 

&>l(X,Y,t) — = " 



2X cos(f cos %t)+(l-X 2 -Y 2 ) sin(| cos ft) 



O (Y VA - ~ 2X Sin( ^ sin 1 *)+( 1 ~- y2 - y2 ) cos (f ft) i 



The expressions (TTOT) . (TTTT) look a bit unwieldy, but this is, apparently, the 
first explicit example of a set of differential equations with smooth coefficients, 
which has attractor of Plykin type on a plane in the Poincare cross-section. 

Figure 3 shows variables X and Y in dependence on time as obtained from 
numerical integration of the differential equations (fTOl) at e=0.72 and K=1.9, 
after exclusion of transients. Visually they look like realizations of random 
processes, as it should be for dynamics on the chaotic attractor. 

Figure 4 shows portrait of the attractor in the three-dimensional extended 
phase space (X, Y, t). To make visible the inherent structure, the picture is 
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presented in the gray-scale technique. Brighter tones correspond to relatively 
larger probability of visiting pixels by orbits on the attractor. Time interval 
on the vertical axis corresponds just to a period of variation of coefficients in 
the equations. In the cross-section of the attractor with the horizontal plane, 
one can observe an object with fractal structure, remarkably similar to that 
discussed in the context of the sphere map (see Fig. 2a). At a qualitative level, 
it may be regarded as an argument in favor of persistence of the hyperbolic 
attractor under modification of the model we undertake. 



X 




Figure 3. Variables X and Y versus time as obtained from numerical integration of 
the differential equations (fTOl) at £=0.72 and if =1.9. The plot relates to sustained 
chaotic regime associated with motion on the attractor; the transients are excluded 




Figure 4. Portrait of attractor of the model ([TO]) at e=0.72 and K=1.9 in the 
extended three-dimensional phase space. The gray-scale technique is used: brighter 
tones correspond to relatively larger probability of visiting the pixels by orbits on 
the attractor 
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Computation of the Lyapunov exponents for the model (fTOl) by means of the 
Benettin algorithm [TT] at e=0.72 and K=1.9 yields Ai « 0.221 and A 2 « 
—0.315. It corresponds to the Lyapunov exponents of the Poincare map Ai = 
X\T pa 0.884 and A 2 = A 2 T pa —1.260. Estimate of the attractor dimension in 
the Poincare section by the Kaplan - Yorke formula is Dky ~ 1.70. 



The hyperbolic nature of the attractor was verified by means of graphical 
representation of manifolds in the Poincare section (like in Fig. 2) and by the 
computations based on the cone criterion (in analogy with Ref. [TT]). As follows 
from those results, attractor is hyperbolic in some parameter range around 
£=0.72 and 7^=1.9. (More details will be published elsewhere.) 



To conclude, this work puts into consideration a non-autonomous flow system 
manifesting chaotic dynamics associated with a hyperbolic strange attractor. 
In the stroboscopic Poincare map it is attractor of Plykin type. In fact, I 
present two versions of the model. In the first version the evolution consists of 
successive stages, and the differential equations have piecewise continuous de- 
pendency of coefficients on time. The Poincare map is expressed analytically, 
as a map on a sphere In the second version, the system is modified in such 
way that it is governed by a set of differential equations with smooth coeffi- 
cients. The modification does not alter the hyperbolic nature of the chaotic 
attractor due to structural stability intrinsic to this object. As the system 
has the minimal phase space dimension required for existence of a hyperbolic 
strange attractor, its investigation, including verification of the hyperbolic- 
ity criteria is essentially simpler, in comparison with models suggested earlier 
as examples of attractors of Smale- Williams type [T0|12fl3| . Appearance of 
concrete examples of systems with hyperbolic strange attractors is of evident 
significance both from the point of view of complementation of mathematical 
concepts with concrete and visible context (see e.g. Ref. [33]), and for ex- 
ploiting these concepts in applications. Hyperbolic chaotic systems may be of 
special interest for applications due to structural stability, that means insen- 
sitivity of the generated chaos to variations of parameters, characteristics of 
elements, technical fluctuations etc. 
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